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1 Abstract 

A simple and efficient computational method is pre- 
sented for unstructured surface grid generation. This 
method is built upon an advancing front techmgue 
combined with grid projection. The projection tech- 
nigue is based on a Newton-Raphson method. This 
combined approach has been successfully implemented 
for structured and unstructured grids. In this pa- 
per, the implementation for unstructured grid is dis- 
cussed. 

2 Introduction 

The first step in obtaining a Computational Fluid Dy- 
namics (CFD) solution is the creation of structured 
or unstructured grids. This step typically requires a 
considerable amount of time and effort on the part of 
a designer. The process of transforming an aerody- 
namic configuration from a Computer-Aided Design 
(CAD) model to a CFD surface grid is referred to as 
the surface-grid generation process. This process is 
often a formidable one for complex geometries such as 
realistic aircraft and spacecraft configurations. CAD 
systems typically represent the surfaces of aerody- 
namic vehicles with a set of parametric surfaces such 
as NonUniform Rational B-Splines (NURBS). Then, 
CFD surface grids are generated on these NURBS 
surfaces. 
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The surface grid can be generated either in a pa- 
rameter space or on approximated/simplified NURBS 
surfaces. Generating grids in a parameter space is 
very common in structured grid generation. This ap- 
proach has two serious restrictions. The first restric- 
tion is that the choice of surface parameterization af- 
fects the CFD surface grid. As shown in ref[l], a poor 
parameterization may cause the CFD surface-grid to 
be highly skewed. There are several ways to allevi- 
ate this problem which have been discussed in great 
detail in [1], The second limitation is that a CFD 
surface grid can not be generated over several over- 
lapping NURBS surfaces. This is the most serious 
restriction. 

In the second method , the NURBS surfaces are 
approximated by a few smaller bi-linear patches. Then, 
the surface grids are generated on these bi-linear patches. 
This method is quite easy to implement, and it avoids 
the problems associated with surface parameteriza- 
tion. However, the resulting grids are close but not 
on the original NURBS surfaces. This problem can 
be alleviated by projecting the resulting grids onto 
the original NURBS surfaces. After the first projec- 
tion, the projected grids may need to be smoothed 
and projected again. The projection techniques de- 
scribed here can be used for structured grid genera- 
tion as well. 

In this study, the NURBS surfaces are first ap- 
proximated by a set of smaller bilinear patches. Then, 
an advancing front technique is used to generate grids 
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on these patches. And, finally these grids are pro- 
jected back onto the NURBS surfaces. In the follow- 
ing sections, the techniques for projecting a point on 
a curve and on a surface are described. Then, a vari- 
ation of the advancing front technique is presented. 
Finally results are summarized. 

3 Projecting on a Curve 



Curves are generally expressed by parametric splines 
such as NURBS as, 

R(u) = {*(«), y(u), r(u)} T , u E [a, 6]. 

Points R(a) and R(b) are the beginning and the end- 
ing points of the curve, respectively. The variable u 
is a parameter which has no geometric significance. 
However, as u increases, the point R(u) always moves 
from the beginning to the end of the curve (Fig. 1). 
This property is referred to as monotone parameter- 
ization which is essential for curve and surface repre- 
sentations. 

The process of projecting a point, r, on a curve, 
R(u), can be achieved by finding a parameter u such 
that the distance between the point, r, and the curve, 
R(u), is minimum and u is constrained to u E [a, 6]. 
This is a constrained minimization problem which 
can be written in terms of distance, d, as 


Figure 1: Parametric Curve 



Figure 2: Projection on a Straight Line 


Two test cases have been selected for this section: 
projecting a point on a straight line and on a NURBS 
curve. 

In the first test case, a point is projected on a 
straight line as shown in Fig. 2. The line is expressed 
in terms of parameter u, as 


d 2 (u.) = f(u.) = | R(u) - r| 2 . (1) 


R(u.) = (R b - R a )u. + R a . (3) 


For curves, the constrained minimization problem is 
relatively simple to solve. A Newton- Raphson method 
is used in this study. In order to find the minimum 
distance, Eq. 1 must be differentiated with respect to 
u and set equal to zero, as 


9f(u) 

du 


3 R(u) 
du 


{R(u) — ?•} = 0. 


( 2 ) 


The projected point, A = R(u), can be determined 
by combining Eqs. 2-3, 

u _ i r ~ Rg) ' {Rb ~ Ra) 

" (R b - R a ) ■ (Rb - Rg)' 

In order to insure that the projected point is on the 
straight line segment, parameter u must be E [0,1]. 
Therefore, the parameter u (calculated from Eq. 4) 
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must clipped as 


u = min{w, 1}, u = max{«, 0}. (5) 


The second test case is based on projecting on a 
NURBS curve. A NURBS curve can bA expressed in 
terms of parameter u as, 


R ( u ) 


E”=o Ni tP (u)wjPi 

E"=o N i, P (u)wi 


( 6 ) 


where Pi are control points (forming a control poly- 
gon), u>i are the weights, and Ni )P are the p-t.h degree 
B-spline basis function dehned on the non-periodic 
and nonuniform knot vector [2]. Combing Eqs. 2 
and 6 yields a nonlinear equation which is solved by 
a Newton-Raphson method. Once the parameter u 
is determined, it has be be clipped as described by 
Eq. 5. This method requires an initial guess which 
can be obtained by sampling the curve, R{u), at sev- 
eral locations. This approach is very robust and effi- 
cient. It takes an average of five iterations to reduce 
the residual by five orders of magnitude. The projec- 
tion process is shown in Fig. 3. Figure 4 shows tin 
results of projecting points onto a NURBS curve. Tin 
dashed-line shows the control polygon for the NURB‘ 
curve. Weights for the first and last control points an 
unity, and weights for the middle control points an 
ten. The solid line shows the NURBS curve, the tri- 
angle symbols show the points before the projection 
and the square symbols show the points after projec 
t.ion. 



A A 

A 



4 Projecting on a Surface 

Figure 4: Projection on a NURBS Curve 

Surfaces are generally approximated by parametric 
surfaces, e.g. NURBS, as, 


R(u) = {*(«), y( u), ;{ u)} T , 
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u = {«i,w 2 } T ,e [(a, 6), (c, rf)], 

where u are the parameters which have no geometric 
significance. However, for a constant u 2 , as ui in- 
creases, the point R(u) always moves from one side 
of the surface to the other side. The process of pro- 
jecting a point, r, on a surface, R( u, ), can be achieved 
by finding u, such that distance, d, between the r 
and R(u) is minimum and u. is constrained to be 
E [(a, b), (c, d,)]. Distance, d, can be written in terms 
parameter u. as, 


d 2 {u) = f(u) = \R{u) - r| 2 . (7) 


In order to find the minimum distance, Eq. 7 must 
be minimized with respect to u. This can be accom- 
plished by setting the gradient of /, V/( t7), equal to 
zero, as 


V/(«) = Gi(S) 


df(u) 

diu 


( 8 ) 


Gi(S) 


dR(u) 

dui 


■ {R(u) 




The above nonlinear system of equations must be 
solved for u. Three test cases are discussed hete: (1) 
projection on a three-dimensional triangle element, 
(2) projection on a bilinear patch, and (3) projection 
on a NURBS surface. 


R, 



Figure 5: Parameter for a Triangular Element 
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Figure 6: Projection on a Triangular Element 


(S2-S 2 )u 1 + (S 1 -S2)u 2 = (S-S 2 ) 


(S 2 -S 1 )u 1 + (S 1 -S 1 )u 2 = (S-S 1 ) 


III + «2 + «3 = 0. 


As shown in Fig. 6, a three dimensional triangle For triangular elements, the above equations are solved 
element can be represented in terms of its parametric for , (7 (see Fig g for definitions of Si , S 2 and S). The 
coordinates as projected point is inside the triangle if 0 < < 1 and 

0 < iii + w 2 + u 3 < 1, otherwise it is outside. If the 
projected point is outside of the triangular element, 
the parameters, must be clipped as described by 
or Eq. 5. 


R(u) = Riui + Rouo + R 3 U 3 , 

111 + «2 + «3 = 1 , 


R(u) — (Ri — + (i? 2 — ^ 3)^2 + R 3 • (9) 


Combining Eqs. 8-9 yields the following system of 
linear equations, 
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The second example is for a bilinear surface, where 
the surface is approximated by a set of structured 
points. A bilinear surface can be decomposed into a 



i+1, j+1 
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Figure 8: Projection Process for a surface 


set of bilinear patches, and each patch, f?(wi,« 2 ), is 
approximated in terms of its parameters (see Fig. 7) 
as 

R(iii , Uo ) = 

( 1 - III )( 1 - U2)Ri,j + ( 1 - Hi )UoRi,j + l + 
«l(l - Un)Ri + l,j + Ulll'jRi + lJ + l- (10) 


Combining Eqs. 8 and 10 yields the following system 
of nonlinear equations, 

— .{R(u)-r} = 0, — .{/?(, J)-r} = 0, (11) 

where, 

OR -> -> -> -> 

~ ' — (1 ll 2 ) ( Ri + 1 j Ri ,j ) T 11 2 ( Ri + 1 J + 1 Ri ,j + 1 ) ■ 

OR -* -* -* -* 

— — (1 ll 1 ) ( Ri J + 1 Ri j ) C ll 1 ( Ri + 1 j' + 1 Ri + 1 ,j ) • 

OH 2 

Equation 11 is solved by Newton-Raphson method. 
It takes an average of 5 iterations to converge. The 
method requires an initial guess which is found by 
sampling the surface at various locations (see Fig. 8). 


The third test case is projecting on a NURBS sur- 
face. A NURBS surface can be expressed as 


R(u) 


E”=0 Ej=0 Ni,r( U l) N 3,q( U 2) W iPi,3 

E "=0 Ej =0 N i,P ( «1 ) N i,q( «2 )ll’i 


Combining Eqs. 8 and 12 yields a system of nonlin- 
ear equations which is solved by the Newton-Raphson 
method. 


5 Advancing Front Technique 

Once the NURBS surfaces are decomposed into smaller 
patches (3/4-sided ), the standard advancing front 
technique is used to generate the interior triangu- 
lation [3]. This triangulation is not on the actual 
NURBS surfaces, but it is close. Then, the result- 
ing triangulation is projected on the actual NURBS 
surfaces as described in the previous sections. The 
VGRID system [3] has been used to generate the ini- 
tial surface triangulation. 

The spacing interpolation in the VGRID system 
is based on a structured background grid [4]. This 
technique simplifies the specifications of grid density 
by introducing nodal and linear sources. The contri- 
butions from nodal sources are inversely proportional 
to the square of the distance. A sample triangulation 
based on a nodal sour®* is shown in Fig. 9. This is 
very similar to the Shepard method of interpolation 
[5]-[6]. However, the contributions from the linear 
sources are modeled similar to the diffusion which 
is not consistent with the nodal sources. The lin- 
ear soured integral-formulation in [4] becomes singu- 
lar near the ends of the elements c Consequently, the 
spacing is more Concentrated near the middle of the 
linear source which results in an asymmetric spacing 
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distribution as shown in Fig. 10-11. In this study, the 
contributions from the linear sources are modified to 
resemble the nodal sources. In order to interpolate 
spacing for a point, P , a point must be found on the 
linear source that is the closest. This can be accom- 
plished by projecting the interpolation point, P , onto 
the linear source as shown in Fig. 12. Equation 4 can 
be used to compute the projected point, a. The as- 
sociated spacing at point a can be calculated as 

S a = (S 2 - S 1 )u + S 1 

where Si and S 2 are the linear source spacings, and u 
is the unidirectional parameter associated with point 
a. For point, P, the contribution from the linear 
source can be treated as a nodal source at location 
at a with spacing S a - Grids based on this formula- 
tion are shown in 13-14. This method can be easily 
extended to curve and surface sources. 

6 Results and Discussions 

Two test cases are projected in this study. The first 
case is a projection onto bilinear surfaces. Figure 15 
shows the results of projecting an unstructured grid 
on bilinear surfaces. Figure. 16 shows the surface 
contours before (solid-line) and after (dash-line) pro- 
jection. Figure 17 shows the projected surface grid for 
an X-15 configuration. This grid has been projected 
onto ten NURBS surfaces. Figure 18 shows the sur- 
face contours before (solid-line) and after (dash-line) 
projection. In both examples the projected grids were 
not distorted. 

Using the projection technique described above, 
it is possible to generate structured and unstructured 
grids on CAD surfaces. The CAD surfaces could be 
overlapping as well. 
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Figure 10: Grid Based on Old Linear Source 
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